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Abstract. - We consider the behaviour of current fluctuations in the one- dimensional par- 
tially asymmetric zero-range process with open boundaries. Significantly, we find that the 
distribution of large current fluctuations does not satisfy the Gallavotti-Cohen symmetry and 
that such a breakdown can generally occur in systems with unbounded state space. We also 
discuss the dependence of the asymptotic current distribution on the initial state of the system. 



Substantial progress in the understanding of nonequilibrium systems has been achieved 
recently through so-called fluctuation theorems [1]. Specifically, the Gallavotti-Cohen fluctu- 
ation theorem (GCFT) can be loosely written as 

P(-(*,t) _ „-at 



p{a, t) 



(1) 



where p(o~, t) is the probability to observe an average value a for the entropy production in 
time interval t and ~ denotes the limiting behaviour for large t. This theorem was first 
derived for deterministic systems [2] (motivated by computer simulations of sheared fluids [3] ) 
and subsequently for stochastic dynamics [4,5]. From [6] onwards there have been successful 
attempts at experimental verification, including for simple random processes such as the driven 
two-level system in [7], Strictly the GCFT is a property of non-equilibrium steady states but, 
for systems with a unique stationary state it is usually also expected to hold for arbitrary initial 
states (see, e.g., [8,9] for discussion on this point). We will refer to this more general property 
of the large deviation function as "GC symmetry" . Some related issues have previously been 
discussed for Langevin dynamics [4,10,11]; we consider the more general case of stochastic 
many-particle systems. 

Specifically, we explore the GC symmetry in the context of a paradigmatic non-equilibrium 
model — the zero-range process [12]. For certain parameter values, this interacting particle sys- 
tem exhibits a condensation phenomenon [13,14] in which a macroscopic proportion of particles 
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pile up on a single site. Condensation transitions are well-known in colloidal and granular sys- 
tems [15] and also occur in a variety of other physical and nonphysical contexts [16]. In [17] it 
was argued that current fluctuations in the asymmetric zero-range process with open bound- 
ary conditions can become spatially-inhomogeneous for large fluctuations — a precursor of the 
condensation which occurs for strong boundary driving. Here, for a specialized case, we ex- 
plicitly calculate the current distribution in this large-fluctuation regime and thus prove a 
breakdown of the symmetry relation (yi . Significantly, we argue that our analytical approach 
predicts that this effect also occurs for more general models. Fianlly, we discuss the relation 
of our results to GCFT breakdowns found in some other works [18-20]. 

Let us begin by defining our model — the partially asymmetric zero-range process (PAZRP) 
on an open one-dimensional lattice of L sites [21]. Each site can contain any integer number 
of particles, the topmost of which hops randomly to a neighbouring site after an exponentially 
distributed waiting time. In the bulk particles move to the right (left) with rate pw n (qw n ) 
where w n depends only on the occupation number n of the departure site. Particles are 
injected onto site 1 (L) with rate a (5) and removed with rate "fw n (f3w n ). If the partition 
function has a finite radius of convergence (i.e, lim n _»oo w n is finite) then for strong boundary 
driving a growing condensate occurs at one or both of the boundary sites [21]. 

We are interested in the probability distribution of integrated current Ji(t), i.e., the net 
number of particle jumps between sites I and Z + 1 in time interval [0, t] . The long-time asymp- 
totic behaviour of this distribution is characterized by the limit of the generating function 

ei(A) = lim -- In (e~ XJl ®). (2) 

t — >oo t 

which implies [5] a large deviation property for the asymptotic probability distribution, pi (j, t) = 
Prob(j; = j, t) , of the observed "average" current ji = Ji/t 

ptV,t) ~ a-**® (3) 

where ei(j) is the Legendre transformation of e;(A), i.e., ei(J) — maxA{ez(A) — Xj}. 

To calculate the current distribution we employ the quantum Hamiltonian formalism [22] 
where the master equation for the probability vector \P t ) resembles a Schrodinger equation 
with Hamiltonian H (see [21] for details). The generating function (e~ XJl ^) can then be 
written as (s\e~ Hlt \Po) where Hi is a modified Hamiltonian in which the terms in H giving 
a unit increase/decrease in Ji are multiplied by e TA [17]. Here \Pq) is the initial probability 
distribution and (s\ is a summation vector giving the average value over all configurations. 
For the current into the system from the left (which can be positive or negative) we consider 
H with lowest eigenvalue §o(A) and corresponding eigenvector |0). In the case where (s\0) 
and (0|Po) are finite, the long-time limiting behaviour is given by 

(e- u » (t) ) ~ (s|0)(6|P )e- SoW * (4) 

In this case we have eo(A) = eo(A) and the form of Hq imposes the GC symmetry relation 

eo(A) = e Q (2E - A) (5) 

which leads, via ®, to the relationship with a — 2Ej and effective field E given by 
e 2E = (af3/~fS)(p/q) L - 1 . The field E can be related to a force F = 2Ek B T. 

While the ground-state eigenvalue calculated in [17] is independent of w n , the latter de- 
termine the form of the eigenvectors (0| and |0). If linv^oo w n is finite then (s|0) diverges for 
some values of A. For a fixed initial particle configuration (0|Pq) is always finite. However, for 
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Table I - Transition values for input current fluctuations in single-site PAZRP 
Values of A Corresponding values of j 

o A l = g i = (P + l-Sf-a-t A = P(j3+-<-6) 2 - al 5 

C — 0+1-6 J» — 13+1-5 ' J b — (/3+ 7 )(/3+7-<5) 

P A 2 _ ()3+7) 2 -Q7-/»+'7 • rj_ 

L — 2 7 « ^ c — ,3 + 7 



3-/33: 2 + A /((5-fl 3: 2 ) 2 +4a7X 2 . _ -(5-fJx 2 ) 



C = 2^ Od 

A 4 _ /3(l-^) + 7 • = Q/37 3 ; 2 -g[/3(l~ 3: ) + 7] 2 . = aj- [/3(1 -z)+7] 2 

L — 7^ - ,e — 2;( / a+7)[/3(l-i)+7j ' Jf — f3(l-x)+-, 



a normalized distribution over initial configurations (e.g., the steady-state) (6|Po) can also di- 
verge (again in the case where w n is bounded) meaning that the asymptotic current distribution 
retains a dependence on the initial state. This has important consequences for measurement 
of the current fluctuations in simulation (or equivalent experiments). Suppose we start from 
a fixed initial particle configuration, e.g., the empty lattice, wait for some time T\ and then 
measure the current over a time interval T^. These are two noncommuting timescales — if we 
take Ti — > oo faster than T± — > oo we will measure the asymptotic distribution of current 
fluctuations corresponding to the fixed initial condition which may differ from the asymptotic 
behaviour of steady-state current fluctuations obtained by taking T\ — * oo before Ti — ► oo. 

We first specialize to the case of the single-site PAZRP, i.e, one site with "input" (left) and 
"output" (right) bonds. In this model explicit calculation of the matrix element (s\e~ Hot \Po) 
is possible. For simplicity we consider here w n — 1, anticipating qualitatively the same effects 
for any bounded w n . We take the case a — 7</3 — 5 in order to ensure a well-defined steady 
state and assume an initial Boltzmann distribution 

oo 

|P > = $> n (l-z)|n) (6) 
n=0 

where \n) denotes the state with site occupied by n particles and the fugacity x = e~^ M < 1. 
The steady state is x = (a + S)/((3 + 7) while x — * gives the empty site. By ergodicity this 
gives the same asymptotic current distribution as any fixed initial particle number. 

Explicit computations yield an integral form for the generating function of input current 

V ' 17 2™ \J Cl {z-\){z-x) 

e -e(z)t X^jux/vx ~ ZVoJifi+j)\ ^ 

c 2 ( z - !)[ z ~ x- 1 ux/vx][z - ux/((3 + 7)] J 

with 

e(z) = a + /3 + j + 5 — vxz — uxz^ 1 . (8) 

Here, for notational brevity we write ux = ae~ x + 6, vx = (3 + je x . and for later use also define 
the parameter combination r) = \/[(f3 + j) 2 — (38 — a"/] 2 — AaffjS. The contour C% (C2) is an 
anti-clockwise circle of radius x + e (e) around the origin of the complex plane with e — > 0. 

In order to extract the large-time behaviour from this integral representation we use a 
saddle-point method, taking careful account of the contributions from residues when the 
saddle-point contour is deformed through poles in the integrand. This yields changes in 
behaviour at the values of A given in Table IJ For 

C ~ 2/303 + 7) {) 
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Fig. 1 - Phase diagram for input current large deviations. Single-site PAZRP with w„ = 1, a = 0.1, 
13 = 0.2, 7 = 0.1, 5 — 0.1. Dotted horizontal line shows mean steady-state current, dashed vertical 
line denotes steady-state initial condition. The symmetry Q is obeyed in the shaded area inside III. 



we find 



eo(A) 



whereas for x > x c we get 



'a(l - e~ A ) +7(1 - e A ) 



a 



a + f3 + j + 5 — 2y/u\v\ 
^a + [3 + j + S — v\x — u\x~ 



A < Ai 
Ai < A < A 2 
A 2 < A < A3 
A 3 < A 



(10) 



e (A) 




u\x~ 



A < Ai 

Ai < A < A 4 

A 4 < A. 



(11) 



We note that the form of eo(A) seen in the regime Ai < A < A3 (A4) is the groundstate 
eigenvalue of Hq [17]. At A = A2 the spectrum of Hq becomes gapless. The changes at Ai and 
A3 (A4) correspond to the divergence of (s|0) and (0|Po) respectively. One immediately sees 
that the symmetry relation JSJ is only obeyed for a limited range of A. 

Via Legendre transformation we obtain the large deviation behaviour of jo — Jo/t. The 
resulting "phase diagram" is shown in Fig. ^ where eo(j) has the following forms in the 
different regions: 



/i(«)7) 
9j 

fl SB. 



( (a-l3-~f+S)(p-5) [3+ 1 - 
\ /3+7-<5 ' a 



7t5 



eo(i) = < ' J V/3+7'/3+7, 



f j (a, 1 )+g j ((3(l-x) + 6(l 

( (l- x ){ a f3 x -5[p{l-x)+ 1 }} 



X 



/3(l-x)+ 7 



I 

II 
III 

IV 

V 
VI 



(12) 
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Fig. 2 - Theory (lines) and simulation (points) for log[p(j,t)/p(—j,t)]. Parameters of Fig with 
initial conditions: (a) x — (empty site), (b) x = 2/3 (steady state). 



with 

fj(a,b) = a + b- Vf+^b + jl j + V] 2 l + Aab (13) 
gj(a,b) = a + jhib. (14) 

The function fj(a,b) is the "random walk" current distribution of a single bond with Pois- 
sonian jumps of rate a to the right and b to the left. The straightline function gj(a,b) gives 
an exponential decay of po(j,t) with increasing j. We now give some brief remarks on the 
physical interpretation of these behaviours. 

In region III, the current across the input bond is dependent on the current across the 
output bond, resulting in a distribution with mean (a[3 — "/8)/((3 + 7) and diffusion constant 
(a/3 + 7<5) /(/3 + 7). In IV (j large and negative) there is a temporary build-up of particles 
on the site (an "instantaneous condensate" [17]) and so to see jo — j, requires a current of j 
across both bonds independently. In I (j large and positive) the piling-up of particles on the 
site means the input bond does not feel the presence of the output bond. The x-dependence 
in region V arises from the possibility of an arbitrarily large initial occupation. II and VI are 
transition regimes involving linear combinations of two different behaviours. They correspond 
to values of A where eo(A) has a discontinuous derivative (cf. a first order phase transition). 
Analogous results for ei(A) and e\(j) which characterize the distribution of outgoing current 
are obtained by the replacements a «-> S, [3 «-* 7, p <-> q, A <-> —A, j <-* — j. 

The GC symmetry states that, e(—j)—e(j) should be a straight line (of slope log[(a/?)/(7<5)] 
in this single-site case) but the results l|12|) imply that this only holds for small j (specifically 
in the shaded region of Fig.^). In Fig. [21 we test this prediction against simulation for both 
input and output bonds. The Monte Carlo simulation results were obtained using an efficient 
event-driven (continuous time) algorithm; for steady-state results the number of histories with 
each initial occupation was weighted according to the known steady-state distribution [21]. 
For increasing measurement times the simulation data converges towards the long-time limits 
predicted by our theory rather than the straight line predicted by GC symmetry. 

Unfortunately, since for increasing times it becomes exponentially more unlikely to measure 
a current fluctuation away from the mean, it is difficult to get long-time simulation data for 
a large range of j. A further check is provided by numerical evaluation of the integral Q) 
followed by numerical Fourier transform to give the finite time distribution of p(j, t) — for small 
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Fig. 3 - Simulation results for log[p(j, t)/p(— j, t] in four-site PAZRP with w n = 1 — 0.5/n, a — 0.1, 
/3 = 0.2, 7 = 0.1, 8 = 0.1, p = 0.525, q = 0.475 and x = 0. Points show data for five bonds at t = 200. 

t this gives excellent agreement with the simulation data; for larger t the integrals converge 
too slowly for the method to be useful. 

We now turn to numerical results for a larger system with a different choice of bounded 
w n , see Fig. |3 In the finite-time simulation regime one again sees indications of violation 
of GC symmetry with bond-dependent form. Physically, we argue that the inhomogeneity of 
the fluctuations across the two different bonds in the single-site PAZRP and the associated 
violation of the GC symmetry is a result of the temporary build-up of particles on the site. 
In general, this possibility is expected to occur in any open-boundary zero-range process 
with liniji^oo w n finite (even when the boundary parameters are chosen so that there is a 
well-defined steady state, i.e., no permanent condensation). 

Mathematically, the observed breakdown of the GC symmetry results from the divergence 
of (s|0) and (0|Po)- For models where the number of particle configurations N is limited, these 
quantities are finite and the relation Q holds for any initial state. However, the limit N — > oo 
does not necessarily commute with the t — > oo limit taken (implicitly) in JQl and (explicitly) 
in J2J. This non-commutation of limits leads in some cases to the violation of QJ even for 
steady-state initial conditions. This and the initial state dependence (due to non-commuting 
timescales) are the main issues highlighted by our work. 

We now give a more general explanation of this GC breakdown and highlight some con- 
nections to previous works. Firstly, consider the observed bond dependence. For systems 
with bounded state space, currents across different bonds differ by finite boundary terms 
which vanish in the long-time limit so any combination of currents has the same large devi- 
ation behaviour. In contrast, for systems with unbounded state space, current fluctuations 
can be spatially inhomogeneous and the boundary terms non-vanishing. However, there is 
always a specific weighted sum of currents for which these boundary terms cancel, giving 
an action functional analogous to heat production (see [5]). For the choice w n = 1 this is 
W = ^^2ii—nEiJi where Ei is the effective field across each bond, e.g., for the single-site 
PAZRP we have e 2E ° = a/7 and e 2El = 8/(3. 

However, it can readily be seen that the large deviations of W still do not satisfy the GC 
symmetry. This is due to the presence of further non-vanishing boundary terms. Consider 
instead the modified action functional (again for w n = 1) 




(15) 
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where {n}(t) represents the configuration of particles at time t and Pq is the initial distribution. 
The fluctuations of this quantity do satisfy the relationship Q even for finite times — this is a 
statement of the transient fluctuation theorem of Evans and Searles [9,23] (see also [24,25] for 
recent experimental tests). Only for bounded state space (finite potentials) do the boundary 
terms containing the initial distribution vanish in the long-time limit leading to recovery of 
the GC symmetry and the steady-state theorem. 

Note that if one measures only a single current (e.g., Jo or J\) but starts with an initial 
distribution corresponding to detailed balance across that bond, the boundary terms cancel 
and the GC symmetry is observed. A particular example is the zero-current case a(3 = "/5 
with an initial equilibrium distribution, x = (a + S) / (f3 + 7) — the current fluctuations across 
both bonds become symmetric e(j) = e(—j) as predicted by the GCFT with E — > 0. This also 
implies the usual Green-Kubo formula and Onsager reciprocity relations [5] . For other values 
of x a breakdown of QJ is still predicted in the E — ► limit (despite the system's ergodicity). 

An analogous apparent breakdown of the GCFT in models with deterministic dynamics 
and unbounded potentials was discussed by Bonnetto et al. [18]. They argue for the restoration 
of the symmetry by removal of the "unphysical" singular terms An earlier study of a model 
with both deterministic and stochastic forces [19, 20] (see [26] for experimental realization) 
found a modified form of heat fluctuation theorem for large fluctuations. In contrast to [18-20] , 
we do not find a constant value for the ratio of probabilities for large forward and backward 
fluctuations. 

A. Rakos acknowledges financial support from the Israel Science Foundation. 
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